A&A manuscript no. 

(will be inserted by hand later] 



Your thesaurus codes are: 
10.05.1;10.06.1;10.07.2;10.08.1 



ASTRONOMY 

AND 
ASTROPHYSICS 



The Old Halo metallicity gradient: 
the trace of a self-enrichment process 



o 
o 
o 

(N 

^■ 
<D • 

ON- 

(N : 



>■ 
i> 

o; 
o. 
o- 

Oh; 



% 



G. Parmentier, E. Jehin, P. Magain, A. Noels and A. A. Thoul* 

Institut d'Astrophysique et de Geophysique, Universite de Liege, 5, Avenue de Cointe, B-4000 Liege, Belgium 

Received date; accepted date 



Abstract. Based on a model of globular cluster self- 
enrichment published in a previous paper, we present an 
explanation for the metallicity gradient observed through- 
out the galactic Old Halo. Our self-enrichment model is 
based on the ability of globular cluster progenitor clouds 
to retain the ejecta of a first generation of Type II Super- 
novae. The key point is that this ability depends on the 
pressure exerted on the progenitor cloud by the surround- 
ing protogalactic medium and therefore on the location of 
the cloud in the protoGalaxy. Since there is no significant 
(if any) metallicity gradient in the whole halo, we also 
present a review in favour of a galactic halo partly build 
via accretions and mergers of satellite systems. Some of 
them bear their own globular clusters and therefore "con- 
taminate" the system of globular clusters formed "in situ" , 
namely within the original potential well of the Galaxy. 
Therefore, the comparison between our self-enrichment 
model and the observational data should be limited to 
the genuine galactic globular clusters, the so-called Old 
Halo group. 

Key words: Galaxy: evolution - Galaxy: formation - 
globular clusters: general - Galaxy: halo 



1. Introduction 

Galactic globular clusters (hereafter GCs) are fossil 
records of the formation of the Galaxy. The understand- 
ing of their formation process would certainly shed light on 
the early galactic evolution. However, at the present time, 
there is no widely accepted theory of GC formation. In 
Parmentier et al. (1999) (hereafter Paper I), we suggest 
a formation scenario based on a self-enrichment process 
such as proposed by Cayrel (1986) and further developed 
by Brown et al. (1991, 1995). 

Our self-enrichment scenario takes place within the 
Fall & Rees (1985) description of the protoGalaxy, namely 
cold clouds embedded in a hot protogalactic background. 
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These cold clouds are assumed to be the progenitors of 
galactic halo GCs. Since they are made up of primordial 
gas, the main advantage of a self-enrichment scenario is 
that it explains in parallel the formation of the clusters 
and the origin of their metal contents. 

The main target of Paper I was to demonstrate con- 
clusively that the gaseous progenitors of galactic halo GCs 
are able to sustain a few hundreds of Type II Supernovae 
(hereafter SNII) without being disrupted. This result is 
in contrast with the widespread idea according to which a 
few supernovae are able to disrupt a Proto-Globular Clus- 
ter Cloud (hereafter PGCC). Furthermore, the large num- 
ber of SNell allowed by our model can explain the amount 
of metals currently observed in galactic halo globular clus- 
ters, and this without any requirement of pre-enrichment 
of the gas. 

The aim of the present paper is to explore further an 
interesting consequence of Paper I, which is also the main 
difference existing between our self-enrichment model and 
the one developed by Brown et al. (1995). The metallic- 
ity that a PGCC can reach through self-enrichment de- 
pends on the pressure exerted by the medium surrounding 
the progenitor cloud and, therefore, on the cloud location 
in the protoGalaxy. The deepest in the protoGalaxy the 
PGCC is located, the highest the final metallicity induced 
by self-enrichment will be. Therefore, we expect to find a 
metallicity gradient throughout the galactic halo. 

The paper is organised as follows. 
In Sect. 2, we briefly review the self-enrichment model pre- 
sented in Paper I, focusing on the link between the final 
metallicity of the PGCC and the pressure exerted on it by 
the surrounding hot protogalactic background. In Sect. 3, 
we examine the different arguments suggesting the exis- 
tence of halo substructures, in order to isolate to which 
one the self-enrichment model can be safely compared. 
In Sect. 4, we compare the model with the observations. 
Sect. 5 explores the putative link between GC metallici- 
ties and their perigalactic distances. Finally, we present 
our conclusions in Sect. 6. 
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2. Self-enrichment model 

According to Fall & Rees (1985), galactic halo GCs were 
formed during the collapse of the protoGalaxy. During this 
collapse, a thermal instability triggers the development of 
a two-phase structure, namely cold and dense clouds in 
pressure equilibrium with a hot and diffuse protogalactic 
background. The temperature of the clouds is assumed to 
remain at a value of 10* K, where the cooling rate drops 
sharply in a metal-free medium. Their masses scale as 
the Jeans mass of a pressure-truncated spherical cloud 
with a temperature ~ 10* K. Since this is of the same 
order of magnitude (~ 10^ M0) as the GC masses. Fall & 
Rees (1985) identify these cold clouds with the progenitor 
clouds of GCs (however, this temperature, and therefore, 
the characteristic mass, is preserved only if there is a UV 
flux able to prevent any H2 formation, the main coolant 
in a metal-free gas below 10* K). As already mentioned, 
the PGCCs are assumed to be metal-free and the forma- 
tion process must therefore explain how the metals are 
provided within each cloud. Within this context, the self- 
enrichment hypothesis was proposed by Cayrel (1986) and 
further developed by Brown et al. (1991, 1995). A first gen- 
eration of stars is assumed to form in the central regions 
of the progenitor cloud. When the massive stars of this 
first generation explode as SNell, all the cloud material is 
progressively swept in an expanding supershell. This su- 
pershell gets chemically enriched with the metals released 
by the exploding massive stars. Since it is a compressed 
layer of gas, it constitutes a dense medium where the for- 
mation of a second generation of stars is triggered. Under 
favourable conditions (see Brown et al., 1995) these sec- 
ond generation stars, formed in the chemically enriched 
supershell, can recoUapse and form a GC. Therefore, the 
first generation SNell provide the GC metals and trigger 
the formation of the GC stars. 

Supernova energetics has been a major criticism of 
the GC self-enrichment hypothesis. However, the main 
target of Paper I was to conclusively demonstrate that 
a PGCC is not necessarily disrupted by SNII explosions. 
Our self-enrichment model is detailed in Paper I. Suffice 
it to say that we compute, for a given hot protogalactic 
background pressure, the supershell velocity during the 
sweeping of the PGCC as a function of the explosion rate. 
Based on this result, we compare the kinetic energy of the 
supershell with the binding energy of the PGCC in order 
to get the maximum number of SNell the cloud is able to 
sustain. 

Whatever the value of Ph , the pressure exerted by the 
hot protogalactic background on the cloud, we find that 
a PGCC can sustain about 200 SNell. Such a large num- 
ber of SNell can provide the amount of metals observed 
in galactic halo GCs. The results of Paper I are summa- 
rized in Table [^ (see Sect. 4.2 for a justification of the 
Ph values). Clearly, GC halo metallicities can be reached 
through self-enrichment. 



Table 1. PGCC masses and metallicities for different val- 
ues of the pressure of the medium confining the PGCCs. 
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Fig. 1. The whole galactic Halo exhibits no clear-cut de- 
pendence of the metallicity on the galactocentric distance. 
It includes different subpopulations, such as the Old Halo 
and the Younger Halo. For the meaning of the different 
symbols, see Sect. 4.1. Data are based on Harris 1996 

We see from Table [^ that the higher the external pres- 
sure is, the higher the metallicity will be. Indeed, since 
the dynamical constraint leads to a constant SNII num- 
ber, assuming a given Initial Mass Function (namely a 
Salpeter one) and a given stellar mass range, a constant 
amount of metals (independent of Ph) is released by the 
first generation massive stars. Since the PGCC mass de- 
creases with increasing external pressure (the Jeans mass 
scales as P^ , Eq. (5) in Paper I), the PGCCs embed- 
ded in a higher pressure medium, namely located deeper 
in the protoGalaxy, reach higher final metallicities. This 
self-enrichment model, contrary to the one developed by 
Brown et al. (1995), implies a metallicity gradient through- 
out the galactic halo. 

At first sight, there is no confirmation of this self- 
enrichment model by the observational data. The galactic 
halo exhibits no significant metallicity gradient (see Fig. y 
data are based on Harris 1996 Q) However, according to 
Zinn (1993), the galactic halo is composed of two distinct 
subpopulations of GCs , what he calls an Old Halo and a 
Younger Halo. The next section presents a review of the 
evidence supporting Zinn's idea. 



^ Data updated 1999 June 22 at http://physun. physic; 
master. ca/ Globular. html. 



3. Two populations of galactic halo GCs 

According to Zinn (1985), the galactic Globular Cluster 
System (hereafter GCS) includes two subsystems: the disk 
GCs I ([Fe/H] > -0.8) and the halo GCs ([Fe/H] < -0.8). 
Both groups also differ in their mean rotational velocity 
and spatial distribution about the galactic center. 
Globular cluster classification has now gone a step further, 
and the halo subsystem could be itself divided into two 
groups. 

3.1. Horizontal Branch morphology 

Zinn (1993) sorts the galactic halo GCs into two groups 
according to their Horizontal Branch (hereafter HB) 
morphology. The HB morphology, namely the colour 
distribution of the stars located along the HB, can be 
described by the index C = {B - R)/{B + V + R) 
where B, V and R are respectively the numbers of blue, 
variables and red HB stars. This index therefore ranges 
from —1 for a purely red HB to +1 for a purely blue 
one. The morphology of the HB is essentially driven by 
the metallicity of the cluster. As such, [Fe/H] is the first 
parameter that governs the HB morphology. However, 
it has long been known that some clusters with similar 
[Fe/H] values present very different HB morphologies 
(e.g. M13 and M3, NGC288 and NGC362). This is the 
so-called second parameter effect: a second parameter 
(at least), in addition to metal abundance, is needed to 
explain the HB morphology. 

Zinn (1993) notices that, while the inner halo GCs 
{D < 8kpc, where D is the galactocentric distance) ex- 
hibit a tight relationship between HB morphology and 
[Fe/H] , the outer halo GCs (D > 8 kpc) show a large scat- 
ter in the same relation. Hence, the requirement of a sec- 
ond parameter to explain the HB morphology is mostly 
needed for GCs located beyond the solar circle. Since, at 
a given metallicity, these clusters have also redder HBs 
than the inner halo GCs, Zinn (1993) divides the galactic 
halo GCs into two groups: the inner halo GCs are classi- 
fied as Blue Horizontal Branch (BHB) clusters, while the 
clusters whose HB types are 0.4 redder than their inner 
halo counterparts (equal metallicity) are labelled as Red 
Horizontal Branch (RHB) clusters. As mentioned above, 
these RHB clusters are mostly located outside the solar 
circle. This link between the infiuence of the second pa- 
rameter and the galactocentric distance indicates that the 
second parameter problem must be related in some way 
to the formation of the Galaxy (Searle & Zinn 1978). 

According to Lee et al. (1993), the most promising 



Due to their high concentration near the galactic center, 
these metal-rich clusters, at least those located within the inner 
~4kpc of the Galaxy, are now thought to be associated with 
the galactic bulge rather than with the disk (Minnitt 1995, 
Cote 1999). 



candidate for this second parameter is age. //age is indeed 
the second parameter, then, at a given metallicity, a RHB 
GC is younger than its BHB counterpart. Therefore, Zinn 
(1993) labels the BHB group and the RHB group respec- 
tively "Old" Halo and "Younger" Halo. To explain these 
putative age differences, and since the "Younger" Halo 
GCs are mostly located in the outer part of the halo, Zinn 
(1993) suggests that the two groups were formed in two 
different ways. The "Old" Halo (hereafter OH) GCs would 
have been formed during the rapid collapse of the proto- 
galactic cloud, while the "Younger" Halo (hereafter YH) 
GCs would have been formed in satellite systems, such 
as dwarf galaxies, which evaded the protogalactic collapse 
and which were later accreted by the Galaxy. 

3.2. GC Ages 

The hypothesis that RHB clusters are younger than their 
BHB counterparts relies on the implicit assumption that 
age is the only second parameter driving the HB morphol- 
ogy, in addition to metallicity. Therefore, it is interesting 
to wonder if age determinations based on the luminosity 
or the colour of the Main Sequence Turn-Off confirm such 
an hypothesis. 

The nature of the second parameter remains a much 
debated question (for contrasting points of view, see the 
reviews by Stetson et al. 1996 and Sarajedini et al. 1997). 
Chaboyer et al. (1996) find that RHB clusters are on av- 
erage 2-3 Gyr younger than BHB clusters, which is consis- 
tent with age being the dominant second parameter. 
This point of view has been reinforced by other stud- 
ies where some clusters have been found to be signifi- 
cantly younger than the bulk of the other galactic halo 
GCs. These are Pal 12 (Stetson et al. 1989), Rup 106 (Buo- 
nanno et al. 1993), Arp2 (Buonanno et al. 1995), 1C4499 
(Ferraro et al. 1995), Pyxis (Sarajedini & Geisler 1996), 
Pal 14 (Sarajedini 1997), Pal5 (Buonanno et al. 1998), the 
outer halo GCs Pal 3, Pal 4, Eridanus (Stetson et al. 1999), 
NCC362, 1261, 1851, 2808 (Rosenberg et al. 1999). All 
these younger GCs are RHB clusters and are located out- 
side the solar circle, where the second parameter effect is 
dominant. Zinn's idea is therefore strengthened. However, 
according to Rosenberg et al. (1999), three RHB clusters 
(NGC3201, NGC5272 and possibly NGC4590) have not 
been proven to be younger. As a result, age appears to 
be one of the most appealing second parameter candi- 
dates (Ferarro et al. 1995, Chaboyer et al. 1996, Stetson et 
al. 1999), but an additional parameter is probably at work 
in a minority of clusters. 



3.3. Kinematic Differences 

Whatever the second parameter is, the concept of two gen- 
uinely distinct halo subsystems gets some support from 
the presence of a kinematic difference between the YH and 
the OH groups. While the mean galactic rotation of the 
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OH group is clearly prograde, it is close to zero, perhaps 
even sligthly retrograde, for the YH group (Zinn 1993, Da 
Costa & Armandroff 1995, Odenkirchen et al. 1997). This 
result has been recently confirmed by the work performed 
by Dinescu et al. (1999). Based on the most complete com- 
pilation of GC proper motions existing up to now, they 
compute the orbits of 38 GCs. Although a sharp trun- 
cation between the different orbital parameters does not 
appear, they show that, on the average, the "Old" Halo 
and the "Younger" Halo GCs exhibit some differences in 
their kinematics and their orbit shapes. The "Younger" 
Halo (RHB) group presents, on the average, a smaller ro- 
tational velocity, larger velocity dispersions, higher orbital 
energies, higher apogalactic distances (-Da > lOkpc) and 
higher excentricities than the "Old" Halo (BHB) group. If 
the RHB group represents an accreted component of our 
Galaxy, the mean rotational velocity suggests that a sig- 
nificant fraction of the outer GCs came from one or more 
ancestral objects on retrograde orbits. 

3.4- Dwarf galaxies accretion 

Dwarf irregular and dwarf spheroidal galaxies, at least 
the most massive ones, also host their own GCS (e.g. 
the Magellanic Clouds, Fornax). Interestingly, in a plot 
of [Fe/H] vs C, the HB morphology index, some of the 
Large Magellanic Cloud (hereafter LMC) GCs fall among 
the outer halo GCs (Da Costa 1993). //age is accepted 
as the dominant second parameter, then these clusters 
are younger than the inner galactic halo GCs. This is 
in agreement with the presence in the LMC of young 
clusters whose masses are within the galactic GC mass 
range (Elson & FaU 1988, Meylan & Heggie 1997). 
The hypothesis that the halo was partly built via 
accretion was underlined by several authors. Lin & 
Richer (1992) suggest that Rup 106 and Pal 12, two 
RHB clusters known to be younger than other GCs with 
similar metallicity (see Sect. 3.2), were tidally captured 
by the Galaxy from the Magellanic Clouds during their 
recent perigalacticon passage. A similar argument holds 
for Pyxis, another RHB and young cluster (Irwin et 
al. 1995). The association between some RHB/younger 
GCs and streams (alignments along great circles over 
the sky which could arise from the disruption of MW 
satellites) is advocated in Majewski (1994) and Fusi Pecci 
et al. (1995). 

All these GCs could therefore have been born well 
apart from the original protoGalaxy, joining our Galaxy 
through later infall events. As such, they are not indica- 
tive of the early formation of the galactic halo. 

Furthermore, Nature currently provides us with an 
example of satellite accretion. The Sagittarius dwarf 
spheroidal galaxy, the closest satellite of the Galaxy, is cur- 
rently undergoing strong tidal distorsions indicating that 
it will probably be disrupted and absorbed by the Milky 



Way (Ibata et al. 1997, Johnston et al. 1999). Based on 
both positional and kinematic data, 4 GCs (M54, Arp2, 
Ter8 with halo metallicities and Ter7 with disk metal- 
licity) unambigously belong to the Sgr dwarf (Ibata et 
al. 1997). These GCs are therefore being incorporated into 
the galactic halo and as such constitute a source of "con- 
tamination" of the genuine galactic GCs, the real tracers 
of the early evolution of the Galaxy. 

Ibata et al. (1997) also strongly suspect the pres- 
ence of a dark halo around the Sgr dwarf. Indeed, dwarf 
spheroidals are among the most dark-matter dominated 
systems known (Mateo 1998) and dwarf irregulars are sig- 
nificanly more dark-matter dominated than are large spi- 
rals (Carignan et al. 1990). The presence of dark matter 
halos around the Milky Way satellites, denser than what is 
found around large spirals, could induce some differences 
in the star and GC formation mechanisms compared to 
what occurs in the protoGalaxy (Larson 1993). 

3.5. Spatial distribution 

Finally, considering the GCs with [Fe/H] < — 1 (in order to 
remove the obvious disk clusters), Hartwick (1987) notes 
that their spatial distribution can be described in terms 
of two subsystems: an inner fiattened distribution and an 
outer more spherical distribution. 

3.6. What does this all mean ? 

The convergence of all the differences mentioned above 
(HB morphologies, ages, kinematic data, galactocentric 
distances and spatial distributions around the galactic 
center) between BHB and RHB GCs adds weight to the 
claim that they form two genuinely distinct groups. The 
existence of two main substructures in the galactic halo 
implies that a hybrid picture could conveniently describe 
its formation (Stetson et al. 1996, Sarajedini et al. 1997, 
Rosenberg et al. 1999). The inner part, populated by 
BHB GCs, would have been formed over a relatively short 
period of time during the collapse of the protoGalaxy 
(Eggen, Lynden-Bell & Sandage 1962), while the outer 
part, which includes most of the RHB GCs, was mainly 
built via accretion and mergers of satellite systems in a 
still ongoing process (Searle & Zinn 1978). In this case, 
the outer halo objects would actually bear little direct 
relevance to the formation history of the main part of 
the Galaxy. As such, they should not be considered when 
comparing our self-enrichment model to the observational 
situation (see Sect. 4.1). 

To disentangle the genuine galactic GCs from those 
formed in satellite systems and accreted afterwards, Zinn 
(1993) suggests to rely on a HB morphology criterion. 
As illustrated above, this approach is indeed fruitful 
since many RHB clusters exhibit peculiarities, such as 
lower ages than their BHB counterparts. However, the 
GCSs of dwarf galaxies are not exclusively composed 



of young GCs and, consequently, the actual situation is 
certainly more complicated. A LMC cluster, Hodge 11, is 
as old as the inner halo GCs (Mighell et al. 1996). The 
Sagittarius cluster system will contribute to both the YH 
(Ter 7 and Arp 2) and OH (Ter 8 and M 54) groups (Da 
Costa & Armandroff 1995). Therefore, the OH subsystem 
may also contain some accreted objects and is not a 
pure sample of GCs formed during the collapse of the 
protoGalaxy main body. NGC2419 and M5 might be 
some of these interlopers. NGC 2419 is the only BHB 
GC located beyond the Magellanic Clouds. While it has 
the same age as M92, an inner halo GC with similar 
metallicity (Harris et al. 1997), it is quite difficult to 
imagine that this metal-poor cluster formed in the inner 
halo and then migrated into the far outer one. The same 
problem stands for NGC 5904 (M5), an outer halo BHB 
GC currently visiting the inner regions of the galactic 
halo (Dinescu et al. 1999). Clearly, Zinn's classification 
needs refinements but nevertheless constitutes a first step 
in understanding the different processes at work during 
the whole halo history. Therefore, in what follows, we 
mainly rely on this BHB/RHB division. 

To close this section, we note that some of the trends 
presented by GCs are also reproduced by field halo stars. 

1. According to Marquez & Schuster (1994), field halo 
stars whose apogalactocentric distances are larger than 
~ lOkpc form a younger group with a larger age dis- 
persion than the inner part of the stellar halo. 

2. Majewski (1992) and Carney (1996) divide the stellar 
halo into a "high halo" sample of stars {{\ZMax\) > 
5 kpc, where Zmux is the maximum height reached by 
the stars above the galactic plane) and a "low halo" 
{{\ZMax\) < 5 kpc). The "high halo", which should be 
dominated by the accreted population, is found in net 
retrograde rotation. In contrast, the "low halo" is in 
prograde rotation. 

3. Finally, Chiba (2000) notes that the stellar halo also in- 
cludes subcomponents characterized by different den- 
sity distributions: the outer part of the halo {D > 
15 kpc) appears to be nearly spherical, whereas the 
inner part exhibits a flattened distribution. 

One interesting point is that the changes in these features 
occur more or less at the galactocentric distance where 
the fraction of RHB clusters increases significantly. The 
stellar halo pattern is therefore consistent with the halo 
GC dichotomy. This gives some support to the hypothe- 
sis of Jehin et al. (1999). Following them, field halo stars 
were, at least partly, formed in GCs from which they es- 
caped through various dynamical processes (disruption or 
evaporation). 



4. The metallicity gradient 



4.I. Importance of the OH/YH division 

The whole galactic GCS exhibits a metallicity gradient 
interpreted as a disk-halo dichotomy, namely, the gradi- 
ent is mainly driven by the high metallicity clusters (disk 
component: [Fe/H] > — 0.8) located within ~8kpc from 
the galactic center (Djorgovski & Meylan 1994). The halo 
group itself ([Fc/H] < — 0.8) presents no clear metallicity 
gradient (see Fig. |l|). However, Sect. 3 provides several ar- 
guments supporting a further meaningful division, mainly 
based on a BHB/RHB (OH/YH) classification, of the halo 
system. As a result, the situation must be reconsidered. 
Following Zinn (1993), the OH GCs are formed during the 
monolithic collapse of the protogalactic cloud while the 
more remote YH GCs are formed in fragments that es- 
cape the protogalactic collapse. Later on, these fragments 
evolve into satellite systems, e.g. dwarf galaxies, bearing 
their own GCs into the galactic halo once they are ac- 
creted by the Milky Way. These dwarf galaxies GCs are 
therefore added to the original galactic GCs. 
According to the model exposed in Paper I, the metallicity 
[Fe/H], induced by the self-enrichment process, is related 
to the hot protogalactic background pressure Ph by (see 
Table 0): 



[Fe/H] = 3.3 + OMogPh 



(1) 



Therefore, to a given pressure distribution of the hot pro- 
togalactic background, Ph{D), corresponds a radial metal- 
licity profile, [Fe/H](Z?). GCs whose progenitor clouds 
were not embedded in this pressure profile will corrupt 
the metallicity gradient if they are considered. 
Since dwarf galaxy GCs were not born in the protoGalaxy, 
their progenitor clouds were not embedded in the same 
pressure profile as the galactic ones. Consequently, they 
must be rejected from any comparison between the theo- 
retical results and the observational data. 

In Fig. |l|, the different types of GCs are marked by 
distinct symbols. OH and YH clusters are respectively rep- 
resented by open and full circles (lists of OH and YH 
GCs are provided in Lee et al. 1994 and Da Costa & 
Armandroff 1995). The crosses label the 3 Sgr GCs with 
[Fe/H] < — 0.8. The open squares ("noHBR" group) stand 
for the halo GCs for which the HB morphology index is 
not given in Harris (1996), probably because no color- 
magnitude diagrams precise enough to define the HB mor- 
phology index were available for these GCs. Indeed, most 
of them are located in the vicinity of the galactic center 
and, therefore, there is no reason to think that they belong 
to the accreted component of the Galaxy. The full triangle 
represents uj Gen, the most peculiar galactic GC (e.g. large 
mass, iron-peak element inhomogeneities in sharp contrast 
with the monometallicity observed in other GCs). These 
peculiarities could be the result of the merger of two GCs 
( Jurcsik 1998) or of the accretion of a dwarf galaxy tidally 
stripped of its stellar envelope (Majewski 2000). 
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Since the self-enrichment model appUes to GCs 
formed in the proto-Milky Way all the progenitor clouds 
are embedded in a common pressure profile and we only 
consider the OH and noHBR groups. This subdivision of 
the halo GCs is important from the point of view of the ex- 
istence or not of a metallicity gradient. In Fig. 0, we show 
[Fe/H] vs logD, limited to the OH and noHBR groups. 
By comparing Figs. II] and 0, we see that the GC metal- 
licity appears to be more strongly correlated with \ogD 
once the presumed accreted component is removed. For 
the whole halo, the linear Pearson correlation coefficient 
is —0.3, corresponding to a probability of correlation of 
99,85%. Considering only the GCs shown in Fig. ||, this 
same coefficient improves to —0.49, corresponding to a 
probability of correlation larger than 99,999%. 
Median values of \ogD in four metallicity bins are given 
in Table 2 for the OH+noHBR groups (the last column 
gives the number of clusters in each bin). The data in this 
table show a monotonic increase of (logD) with decreasing 
metahicity over the range [Fe/H] ~ -0.8 to [Fe/H] ~ -2.4. 



Table 2. Metallicity vs logD for the OH-^noHBR groups. 



[Fe/H] 




(logD) 


n 


-1.2 < [Fe/H] < - 
-1.6 < [Fe/H] < - 
-2.0 < [Fe/H] < - 
-2.4 < [Fe/H] < - 


-0.8 
-1.2 
-1.6 
-2.0 


0.41 
0.64 
0.89 
1.02 


13 

28 
25 
9 



On the average, the more remote GCs are more metal- 
poor. This is in agreement with the self-enrichment model 
where the metallicity gradient is due to the decrease of the 
pressure exerted by the hot protogalactic background on 
the PGCCs as the galactocentric distance increases (see 
Table H). To compare the theoretical metallicity gradient 
to the observational situation, one still needs an expression 
for the pressure profile Ph{.D) (see Eq. (|l|)). 

4-. 2. Pressure profile of the hot protogalatic background 

Concerning this point, the situation is rather complex 
since there is obviously no agreement about the scaling 
of the Ph{D) relation in the literature. 

The luminous components of galaxies form through 
the collapse of gas in gravitationally dominant halos of 
dark matter. The density profile of such a halo is conve- 
niently described as a singular isothermal sphere (White 
& KauflFmann 1994): 



V' 



1 



(3) 



and 

^ o GMID) 

V^ = -^ = const. 

D 

In these equations, Vc is the circular velocity of the gas in 
the dark matter potential well of the Galaxy and G is the 
gravitational constant. Since we are mainly interested in 
the proto-Milky Way, the circular velocity Vc is taken to 
be 220kms-i (FaU & Rees 1985). 

Two timescales are important in determining the fur- 
ther evolution of the gas component within this halo of 
dark matter : 



1. the dynamical (or free-fall) time: 

_ D 

Tdyn — T^, 

2. the radiative cooling timescale: 
3 utkT 



'^cool 



2 rijUpA 



(4) 



(5) 



where n^, n,, and rit are respectively the ionic, electron 
and total number densities, k is the Boltzmann con- 
stant, T is the gas temperature, and A is the cooling 
function (Sutherland & Dopita 1993). 

When Tdyn < Tcooh the gas undergoes quasi-static 
contraction, but once TcooI < Tdyn, the gas cloud can no 
longer maintain itself in quasi-static equilibrium and the 
collapse proceeds on a free-fall timescale (e.g. Rees & 
Ostriker 1977). 

According to Fall & Ress (1985), during the collapse 
of the protoGalaxy a two-phase structure, namely cold 
clouds embedded in the remaining hot protogalactic back- 
ground, grows in the collapsing gas. The hot component 
is depleted by the condensation of the cold clouds until 
its cooling and dynamical timescales are comparable. The 
hot gas is expected to remain near the virial temperature 
of the halo described by Eq. (g) : 

n^^fihniH ^1.7 X 10'' K (6) 

where /i^ is the mean molecular weight (~ 0.6 for a pri- 
mordial ionized plasma, with X ~ 0.76 and Y ~ 0.24). 
At this temperature, the value of the cooling function A 
is 4.8 X 10^24 ergcm^s-i (Sutherland & Dopita 1993). 
Combining Eqs. (0-0)' ^'^ have 



Ph = 5.3 X 10-i"Z?,tJ, dyne cm 



kpc 



(7) 



where Dkpc is the galactocentric distance expressed in kpc. 
The corresponding metallicity gradient is therefore 
(dashed line in Fig. g): 



[Fe/H]- -1.34- 0.5 logi^fcp. 



(8) 



^^^^-^D^ 



(2) 



In the Fall & Rees model, the hot component of the gas 
in which the PGCCs are embedded is at the virial tem- 
perature and has a density such that TcooI — Tdyn- The 



model proposed by Murray & Lin (1992) is slightly differ- 
ent. It relies on the TcooI — Tdyn condition and also on the 
hypothesis of hydrostatic equilibrium for the hot gas: 



1 dPfe 
Ph dD 



y2 

^ c 

D ' 



(9) 



The peculiarity of their model is that the rotation curve 
of the collapsing protoGalaxy is not assumed to be flat, 
but rather slightly decreasing with galactocentric distance. 
They argue that since the central regions, with the highest 
density, are the first to collapse, followed by regions from 
larger initial radii, the gravitational potential may be more 
centrally peaked than today. They adopt Vc{D) ex Z?"^/'*. 
It should be noted that the value of the exponent ensures 
that the temperature of the hot gas obeys to Eq. (0) with 
Vc and Th depending on D (their Eq. (3.9)). In contrast 
with Fall & Rees (1985), their pressure profile for the hot 
protogalactic background scales as D^^: 



Ph = 1.25 X lO-^D^Tp^^ dynccm- 



(10) 



The corresponding metallicity gradient is therefore steeper 
(plain line in Fig. ||) than in Eq. (||): 



[Fe/H] = -1.15-logDfcpc 



(11) 



A third approach is suggested by Harris & Pudritz 
(1994). According to them, it is natural to expect that 
the hot and diffuse phase of the protoGalaxy is at the 
virial temperature and in hydrostatic balance with the 
isothermal dark matter potential well. Under these as- 
sumptions, the pressure distribution scales again as D^^ , 
Tdyn and Tcooi scale respectively as D (Eq. ^) and as D"^ 
(Eq. (P)). There is thus a "critical" galactocentric distance 
Dcrit such that 

- Tcool < Tdyn when D < Dcrit'- 

no quasi-static equilibrium is possible; 



Tcool — Tdy 



when D > Dr 



the gas can remain in hydrostatic equilibrium. 

In contrast to Fall & Rees (1985) and Murray & Lin 
(1992) models, there is no equality between TcooI and 
Tdyn over the whole galactic halo and, therefore, this 
property can no longer be used to determine Ph{D). 
However, it is interesting to note that the steady rise 
of the ratio T^ooi/Tdyn with galactocentric distance is a 
property found in the halos of giant ellipticals, e.g. M 87 
(Fabricant et al. 1980). 

Murray & Lin (1992) and Harris & Pudritz (1994) 
give the same scaling law for the pressure distribution 
{Ph oc D~'^), but not from the same hypotheses! This 
illustrates the difficulty to derive an expression for the 
metallicity gradient from the self-enrichment model. 
In what follows we will use the results given by Eqs. (^ 
and (llT 




Fig. 2. Comparison between 2 theoretical [Fe/H] vs logl? 
relationships and the metallicity radial distribution func- 
tion of the OH-fnoHBR groups. The plain and dashed 
curves represent the self-enrichment model combined with, 
respectively, the Murray & Lin (1992) and the Fall & Rees 
(1985) pressure profiles. GC symbols are identical as in 
Fig. I 



4-3. Discussion 

In Fig. pi we compare the results given by Eq. (0) (dashed 
line) and Eq. (O) (plain line) to the observational data. 
The region of interest, namely where the bulk of the 
OH group is located, is the galactocentric range 1 to 
30kpc. The metallicity intervals predicted by Eqs. (Isl) 
and (11) for these values of the galactocentric distance 
are respectively [— 2dex, — 1.35dex] and [— 2.65dex, 
— 1.15dex]. The second model (self-enrichment model 
combined with Ph oc D~^, in this case the Murray & Lin 
(1992) pressure distribution) is a better description for 
the observed metallicity range (see Table ||) . 

The theoretical metallicity gradients corresponding to 
Eqs. (|) {Ph oc £>-!) and (|n]) {Ph oc D-^) are respec- 
tively Ape/HJ/AlogD = -0.5 dex and -Idex. It is not 
straightforward to compare these results with the obser- 
vational gradient because of the rather high dispersion in 
the [Fe/H] vs \ogD plot, partly due to measurement errors 
in [Fe/H] and logl?. The observed dispersion can also be 
explained by the GC orbital motions. Indeed, our model 
predicts a relation between the GC metallicities and the 
galactocentric distances of their formation site. But the 
GCs were carried away from their formation sites through 
their orbital motions. The initial radial distribution of 
globular cluster abundances has therefore been modified 
but we can only use the current one. A least-squares fit to 
the OH-|-noHBR groups, which takes into account the un- 
certainties in both coordinates (Press et al. 1992), yields: 



[Fe/H] = (-1.35 ± 0.26)logi:» + (-0.64 ± 0.24). 



(12) 
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If the noHBR clusters are not included, the slope of the 
least-squares fit is even steeper with a value of —1.68 ± 
0.41. The errors on [Fe/H] are assumed to be ±0.15 dex 
(King 1999). The errors on logl? are deduced from the 
comparison between the observational points and a "clas- 
sical" least-squares fit to the \ogD vs [Fe/H] plot with 
ciogD=l- This leads to an error of ±0.4 in logD. We cau- 
tion that, by doing so, the errors are assumed to be nor- 
mally distributed. This is probably not the case for logD, 
since the galactocentric distances of the formation sites are 
replaced by the current galactocentric distances. There- 
fore, the dispersion in logD cannot be attributed solely to 
measurement errors. The slope of Eq. (|lj) is consistent 
with the one of Eq. ( [ll|) at the l.Scr level, while it is dif- 
ferent from the slope given by Eq. (g) at the 3.3ct level. 
As a conclusion, Eq. (^ better describes the observed ra- 
dial distribution of GC metallicities when it is combined 
with a pressure profile scahng as D~^ (Eq. (|ll|)) rather 
than D~^ (Eq. (|8|)). However, Eq. ( [ll| ) somewhat under- 
estimates the mean observed metallicity at a given galac- 
tocentric distance. This can be due to: 

1. an underestimate by the self-enrichment model for the 
mass of metals ejected by the SNell in the PGCCs 
(uncertainties in, e.g., the SNII yields, the high- mass 
stellar mass spectrum; see Paper I for a discussion of 
the second point); 

2. an underestimate of the hot protogalactic background 
pressure; 

3. the possibility that GCs formed deeper in the galactic 
potential well than their current location in the halo 
(see Sect. 5 for a related point). 



In addition to the metallicity gradient, the model 
presents a second interesting consequence. According 
to Sarajedini et al. (1997), once the search for an age- 
metallicity relationship is restricted to the Old Halo 
group, there is no evidence for such a relation. This 
result is confirmed by the careful work carried out by 
Rosenberg et al. (1999). Their conclusion is consistent 
with a single mean value for the age, independent of 
[Fe/H]^ Clearly, a parameter other than the age is needed 
to explain the metallicity range observed in the OH/BHB 
group. The self-enrichment model provides us with such 
a parameter, namely the external pressure around the 
PGCCs. Combined with a pressure distribution scahng 
as D~^, the model explains fairly nicely the galactic halo 
GC metallicities, without any requirement for an age- 
metallicity relation (in the usual sense of age decreasing 
with metallicity). The final proportion of heavy elements 
is mostly fixed by the pressure exerted by the medium 
surrounding the progenitor clouds. 



^ Again, the same pattern is observed for the field stars: there 
is no age-metallicity correlation for the FHS (Marquez & Schus- 
ter 1994). 



If, in the future, the pressure distribution was proven to 
scale as D~^ instead of D~^, then this would mean that 
our model underestimates/overestimates the amount of 
chemical enrichment in the inner/outer regions of the halo. 

Finally, the steady decline of the hot gas pressure with 
galactocentric distance should also lead to a gradient in 
the mass of the PGCC (see Table p. That such a mass 
gradient could still be observable today is not a certainty. 
Indeed, the initial mass of the PGCCs undergoes changes 
due to various evolutionary processes: 

a. The formation of the second stellar generation: there 
is no reason why this star formation episod would take 
place with the same star formation efficiency in each 
cloud. 

b. Mass loss (evaporation): once they are formed, GCs 
undergo mass loss due to their interactions with the 
galactic gravitational field (Meylan 2000). 

Evaporation and random star formation efficiency are 
sources of scatter in the relation between the PGCC and 
GC masses. In addition, mass determinations of GCs are 
still uncertain at least by a factor 2 (Meylan 2000). There- 
fore, and as pointed out by Brown (1993): "Even if the 
Jeans mass did play a role in the formation of GCs, it is 
unlikely that a one-to-one relationship exists between the 
proto-cluster cloud mass and the present cluster mass." 
GC absolute visual magnitudes My could also be used 
to search for a relic of the mass gradient since obtaining 
GC integrated luminosities is much easier than determin- 
ing GC masses. However, the absolute visual magnitude is 
not an accurate estimation of the GC mass because of the 
scatter introduced by the mass-to-light ratios in the mass- 
luminosity relation. For instance, in the Pryor & Meylan 
(1993) compilation, the mass-to-light ratios of most of the 
OH GCs range from 1.1 to 3.8. Despite these various un- 
certainties, Brown (1993) finds a steady decline in My for 
D < 15kpc in a plot of average My versus average logD 
for the galactic halo clusters. Therefore, his conclusion fur- 
ther supports the hypothesis of an origin of galactic halo 
GCs in clouds having the Jeans mass. 



5. Were GCs formed near their perigalacticon? 

The initial halo metallicity gradient was severely altered 
by the accretion of GCs formed in dwarf and/or irregular 
galaxies. The genuine Milky Way GCs moving away from 
their formation sites is another cause of alteration. It 
would be very useful to know the galactocentric distances 
at which the GCs were formed, in order to better evaluate 
the validity of Eqs. (|) and ([l|). 

In Sect. 4.3, it was suggested that GCs might have 
formed deeper in the galactic halo than their current 
location. However, we caution that this result relies on 
the self-enrichment model and on the assumed pressure 
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distribution. Therefore, this outcome cumulates uncer- 
tainties from both and it certainly does not stand on a 
firm support. Even so, previous papers already suggest 
that the correlation between [Fe/H] and logZ?p, where 
Dp is the perigalactocentric distance of the GCs, may be 
stronger than between [Fc/H] and logD. 

Freeman & Norris (1981) notice "a clear gradient in 
the [Fe/H]-I?p plane" . According to them, this may be a 
hint that the clusters did form near perigalacticon. Van 
den Bergh (1995) also notes that the GC metallicity cor- 
relates somewhat more strongly with Dp than it does with 
the present GC galactocentric distances. Nevertheless, the 
method used to derive these perigalactic distances is some- 
what problematic. It assumes that the current GC tidal 
radius is mostly set by the galactic tidal field at the closest 
approach of the cluster to the galactic center, namely the 
perigalacticon (King 1962). As a result, tidal radii of well 
observed GCs are used to estimate their perigalactic dis- 
tances. By so doing, the GC internal processes, especially 
the two-body relaxation, are neglected (Meziane & Colin 
1996). These internal processes lead to a replenishment of 
the outer regions of the cluster between two perigalactic 
passages and therefore modify the outer radius set at the 
perigalacticon. Thus, it is dangerous to rely on the cur- 
rent outer radius of GCs to derive their perigalactocen- 
tric distances. In order to avoid this problem, perigalactic 
distances derived from the computation of GC orbits (Di- 
nescu et al. 1999) are used in Figs. || (all GCs with known 
perigalactic distances) andH (OH GCs). Another cluster, 
NGC 6522 (Terndrup et al. 1998), is added to Dinescu's 
list. The sample is by far smaller than in Van den Bergh 
(1995). This is due to the necessity to know the proper 
motions in order to compute the GC orbits in a given 
galactic potential. 

The (logDp, [Fe/H]) plot (Fig. Q) does not appear tighter 
than the (logD, [Fe/H]) one (Fig. ||). The linear Pearson 



correlation coefficient in Fig. ||is only —0.31, correspond- 
ing to a probability of correlation of the order of 90 %. In 
contrast, considering the (logD, [Fe/H]) plot for the same 
sample of GCs, the correlation coefficient is still —0.43, 
corresponding to a probability of correlation of 96 %. The 
orbits only provide evidence that the more metal-rich halo 
clusters ([Fe/H] > —1.4) are concentrated towards the 
galactic center {Dp < 3kpc). Therefore, these new data 
do not really confirm the suggestion made by Freeman & 
Norris (1981). New perigalactic distances would be help- 
ful to give a definitive answer to the existence of a link 
between the perigalactocentric distances and the mctal- 
licities of galactic halo GCs. 

6. Conclusions 

The self-enrichment model of galactic halo GCs (see Pa- 
per I for a detailed description) has been compared to the 
observational situation and the conclusions are as follow: 

1. The final metallicity induced by the self-enrichment 
process depends on Ph, the pressure exerted on the 
PGCC by the hot protogalactic background medium. 
The result is in agreement with galactic halo GC 
metallicities. 

2. There is a range of halo GC metallicities due to the 
variations in Ph with the galactocentric distance. Con- 
sidering the clusters located between 1 and 30kpc, 
the ranges of metallicities [Fe/H] are 0.75 and 1.50 dex 
when the self-enrichment model is combined with Ph oa 
D~^ and Ph oc D~^ respectively. The second solution 
is in better agreement with the observed metallicity 
range. Moreover, when the Murray & Lin (1992) pres- 
sure distribution is used, the theoretical metallicity in- 
terval over the galactocentric range 1 < Dkpc < 30 is 
—2.65 < [Fe/H] < —1.15. This result is in nice agree- 
ment with the observations in the galactic halo GCs. 
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3. Because of the expected decrease in Ph with increas- 
ing galactocentric distance, the model induces a metal- 
licity gradient throughout the galactic halo. Such a 
metallicity gradient is indeed observed once the GCs 
suspected to have been accreted by the Milky Way are 
removed. These clusters were born in fragments which 
evaded the initial protogalactic collapse and have ex- 
perienced their own chemical evolution. The division 
between the galactic and accreted components of the 
halo is mainly based on the BHB/RHB (OH/YH) clas- 
sification introduced by Zinn (1993). Indeed, it seems 
likely that the halo GCs consist of clusters with more 
than one origin. Since the accreted clusters did not take 
part in the formation and early evolution of the galac- 
tic halo, their progenitor clouds did not share the same 
external pressure distribution as the genuine galactic 
proto-GCs. Therefore, they should not be taken into 
account in the self-enrichment model. 

Again, the observed radial distribution of GC abun- 
dances favours a background pressure profile scaling 
as D~^ rather than D^^. 

However, it should be noted that the scatter of the 
data about the model lines in Fig. g (and also in 
Fig. m exceeds the observational uncertainties. As al- 
ready stated, this can be due to the GC orbital motions 
which carry them away from their formation sites. Fur- 
thermore, it is certainly an oversimplification to con- 
sider that Ph is the only parameter determining the 
GC metallicity. Other parameters must interfere (e.g. 
stellar mass range, SNII yields, ...). 

4. In our model, no age-metallicity relation is required to 
explain the different GC metallicities. Actually, there 
is no compelling evidence for an age-metallicity rela- 
tionship among halo GCs (Buonanno et al. 1998), and 
especially once the sample is limited to the BHB/OH 
group (Sarajedini et al. 1997). In this group, all GCs 
are coeval according to Rosenberg et al. (1999). This is 
in agreement with our self-enrichment model where we 
see an enhanced chemical enrichment with decreasing 
galactocentric distance rather with time. 
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